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Magnetic nanostructures on non-magnetic or magnetic substrates have attracted strong attention 
due to the development of new experimental methods with atomic resolution. Motivated by this 
progress we have extended the full-potential Korringa-Kohn-Rostoker (KKR) Green function method 
to treat non-collinear magnetic nanostructures on surfaces. We focus on magnetic 3d impurity 
nanoclusters, sitting as adatoms on or in the first surface layer on Ni(001), and investigate the 
size and orientation of the local moments and moreover the stabilization of non-collinear magnetic 
solutions. While clusters of Fe, Co, Ni atoms are magnetically collinear, non-collinear magnetic 
coupling is expected for Cr and Mn clusters on surfaces of elemental ferromagnets. The origin 
of frustration is the competition of the antiferromagnetic exchange coupling among the Cr or Mn 
atoms with the antiferromagnetic (for Cr) or ferromagnetic (for Mn) exchange coupling between the 
impurities and the substrate. We find that Cr and Mn first-neighbouring dimers and a Mn trimer on 
Ni(001) show non-collinear behavior nearly degenerate with the most stable collinear configuration. 
Increasing the distance between the dimer atoms leads to a collinear behavior, similar to the one 
of the single impurities. Finally, we compare some of the non-collinear ab-initio results to those 
obtained within a classical Heisenberg model, where the exchange constants are fitted to total 
energies of the collinear states; the agreement is surprisingly good. 



I. INTRODUCTION 

Theoretically, extensive work is carried out in the area 
of complex non-collinear magnetism, particularly for sur- 
face and bulk systems. In fact a lot of interesting physics 
would be missed if only collinear magnetic structures 
were considered. In fact, magnetic nanostructures on 
magnetic or nonmagnetic substrates are attractive to 
the scientific community due to their novel and unusual 
propertiesi*2iii4 being of relevance both for theory as well 
as for the applications in the magnetoelectronics devices. 

One of these properties is the non-collinear magnetic 
order occurring for geometrically frustrated antiferro- 
magnets, e.g. on a triangular lattice, in disordered sys- 
tems, exchange bias systems, and molecular magnets, or 
for systems which exhibit either competing exchange in- 
teractions, or between exchange and spin-orbit interac- 
tions. A simple model for frustration is the following: 
Starting with an antiferromagnetic (AF) Cr dimer, the 
addition a third Cr atom to form an equilateral trian- 
gle leads to a frustrated geometry. Each atom would 
like to couple AF to both other atoms. Since this is im- 
possible, the moments of the three atoms rotate until 
a compromise is found. The ground state is then non- 
collinear, characterized by an angle of 120° between each 
two atoms. The same situation will also occur for the AF 
Cr dimer, since the interaction of both Cr atoms with the 
ferromagnetic substrate atoms is either ferromagnetic of 
antiferromagnetic. As we will show in this paper, also in 
this case a non-collinear structure can result. 

The majority of the ab-initio methods available for the 
treatment of non-collinear magnetism make explicit use 
of Bloch's theorem and are thus restricted to periodic 
systems (bulk or films). Then, even for collinear mag- 
netism, one needs huge supercells to simulate impurities 



in a given host (bulk or film) in order to avoid spurious 
interactions of the impurities from adjacent supercells. 
A few methods have been developed to treat free clus- 
ters, but to our knowledge no ab-initio methods exist for 
the investigation of non-collinear magnetism of clusters 
in bulk or deposited on surfaces. 

First non-collinear calculations by the KKR Green 
function method, though not self-consistent, were al- 
ready performed in 1985. Oswald et al~ could show by 
using the method of constraints that the exchange in- 
teraction between the moments of Mn and Fe impurity 
pairs in Cu is in good approximation described by the 
cos ^-dependence of the Heisenberg model. 

Sandratskii et al& and Kiibler et alA& pioneered the 
investigation of non-collinear magnetic structures using 
self-consistent density functional theory. One of the first 
systems studied by Sandratskii et al& was the spin spiral 
of bcc Fe with the KKR method. Later on, A-Fe was a 
hot topic, and the appearance of the experimental work 
of Tsunoda et alJLlii led to the development of other first- 
principles methods able to deal with non-collinear mag- 
netism such as LMTOii, ASWi2 and FLAPW. 13 i 14 i 15 

Several papersi&i! describe how symmetry simplifies 
the calculational effort for the spiral magnetic structures 
in the case of perfect periodic systems — this involves the 
generalized Bloch theorem. In ab-initio methods, this 
principle is used together with the constrained density 
functional iheoryi^i^ giving the opportunity of studying 
arbitrary magnetic configurations where the orientations 
of the local moments are constrained to nonequilibrium 
directions. 

Concerning free clusters, few methods are developed. 
For example, Oda et alm& developed a plane- wave pseu- 
dopotential scheme for non-collinear magnetic structures. 
They applied it to small Fe clusters for which they found 
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non-collinear magnetic structures for Fes and linear- 
shape Fe3 . This last result was in contradiction with the 
work of Hobbs et alm^ who found only a collinear fer- 
romagnetic configuration using a projector augmentcd- 
wave method. Small Cr clusters were found magnetically 
non-collinearj2£ as shown also by Kohl and BertschS^ us- 
ing a relativistic nonlocal pseudopotential method after 
optimization of the ionic structure by a Monte Carlo 
technique. However, within the generalized gradient ap- 
proximation (GGA) of density-functional theory, Hobbs 
et alM> find that in many cases the non-collinear states 
can be metastable, while the ground-state solutions are 
collinear and arise after geometrical optimization of the 
free-standing clusters. 

One main result of Oda et alM- and Hobbs et a/i2I con- 
cerns the variation of the magnetization density with the 
position. The spin direction changes in the interstitial 
region between the atoms where the charge and mag- 
netization densities are small, while the magnetization is 
practically collinear within the atomic spheres. This sup- 
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Here, E is the energy and R n , R n ' refer to the atomic 
nuclei positions. By r< and r> we denote respectively 
the shorter and longer of the vectors r and f* which de- 
fine the position in each Wigner-Seitz cell relative to the 
position R n or R n i . The wavefunctions R'l (f ; E) and 
H2(r;E) are, respectively, the regular and irregular so- 
lutions of the Schrodinger equation for the potential V n at 
site n, being embedded in free space; L = (I, m) is a com- 
bined index for angular momentum quantum numbers; I 
is truncated at a maximum value of /max- The first term 
on the LHS of Equation JQ) is the so-called single site 
scattering term, which describes the behavior of an atom 
n in free space. All multiple- and back-scattering infor- 
mation is contained in the second back- scattering term 
via the structural Green functions G^^E) which are 
obtained by solving the algebraic Dyson equation: 

Gl n L>(E) = °GTl>{E) 

+ y, Gii;,(E)Aq;, L ^E)Gi:,?i,(E) (2) 

n",L"L"' 

Equation J5J follows directly from the usual Dyson eq. of 
the form G = G + G AV G, with AV the perturba- 

o 

tion in the potential and G the reference system Green 
function. The summation in (0) is over all lattice sites 
n" and angular momenta L" for which the perturbation 

At^, L ,„(E) = il„ L ,„{E) - i^, L ,„(E) between the t ma- 
trices of the real and the reference system is significant 



ports the use of a single spin direction for each atomic 
sphere as an approximation in order to accelerate the 
computation; this approximation is followed also here. 

The aim of this work is to present a method based 
on the full-potential KKR scheme 23 which can deal with 
non-collinear magnetism in systems of reduced symme- 
try. This method is ideal for treating impurities or small 
clusters on surfaces or in bulk. As an application we 
study small 3d clusters on the Ni(OOl) surface where we 
find complex magnetic configurations. 



II. NON-COLLINEAR KKR FORMALISM 

The KKR method uses multiple-scattering theory in 
order to determine the one-electron Green function in a 
mixed site and angular-momentum representation. The 
retarded Green function is expanded as: 



(1) 

I 

(the i-matrix gives the scattering amplitude of the atomic 

o 

potential). The quantities G™ l L ,(E) are the structural 
Green functions of the reference system. For the cal- 
culation of a crystal bulk or surface, the reference sys- 
tem can be free space, or, within the tight-binding KKR 
formulation^ a system of periodically arrayed repulsive 
potentials. After the host (bulk or surface) Green func- 
tion is found, it can be used in a second step as a reference 
for the calculation of the Green function of an impurity 
or a cluster of impurities embedded in the host. 

The algebraic Dyson equation © is solved by matrix 
inversion, as we will see later on in Equation (|17l) . In 
case of spin-dependent electronic structure, spin indexes 
enter in the i-matrix, the Green functions and in Eq. J2J. 
Especially in the case of non-collinear magnetism, these 
quantities become 2x2 matrices in spin space, denoted 
by t and G. 

Once the spin-dependent Green function is known, all 
physical properties can be derived from it. In particular, 
the charge density n(f) and spin density m(r) are given 
by an integration of the imaginary part of G up to the 
Fermi level Ep and a trace over spin indexes s (putting 
the Green function in a matrix form in spin space): 

n{r) = --lmTr s J G(r,f;E)dE (3) 
m(r) = ImTr s / a G(r, r; E) dE. (4) 



-iVEj2 R L^<;E)H r l(r>;E)6 nn , + £ R n L (r; (E)R£(f> '; E) 

L LL' 
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Here, <r = ((T x ,<T y ,a z ) are the Pauli matrices and Tr s 
means the trace operation in spin space. 

The basic difference between non-collinear and 
collinear magnetism is the absence of a natural spin quan- 
tization axis common to the whole crystal. The density 
matrix is not anymore diagonal in spin space as in the 
case of collinear magnetism. Instead, in any fixed frame 
of reference it has the form 



p(r) = 



= \ Mr) + d ■ rn{f)] (5) 



At any particular point in space, of course, a local frame 
of reference can be found in which p is diagonal, but this 
local frame can change from point to point. 

In order to deal with non-collinear magnetism, we have 
to solve the appropriate Dyson equation. First we define 
the reference system which is a perfect surface character- 
ized by collinear magnetism. Although the collinearity 
of the reference system is not a necessary requirement, it 
serves our purpose of calculating the electronic structure 
of the ferromagnetic or nonmagnetic surfaces which are 
used as reference systems. Thus the host Green functions 

o o 

G and i-matrices t are assumed diagonal in spin space. 
In this way, in the case of a magnetic host, a global spin 

o o 

frame of reference is defined. The host G and t are thus 
of the form: 



G(E) 



G n (E) 
. G U (E\ 



;t(E) 



t n (E) 
. t u (E)j 



(6) 



Then the perturbed system is constructed. The impurity 
atoms which might couple magnetically in a non-collinear 
way reside on the surface, perturbing the potential at a 
few neighboring sites (atoms or empty cells representing 
the potential in the vacuum). Within this finite cluster 
of perturbed sites the magnetization can be non-collinear 
leading to the appearance of non-diagonal elements of the 
i-matrix: 



t(E) 



ttf(E) t n (E) 
t iT (£) t u (E) 



(7) 



The non-diagonal i-matrix contains the information on 
spin-flip scattering by the atomic potential. 

At this stage an approximation enters our method. It 
is assumed that, separately for each atom, there exists 



an intra-atomic spin quantization axis common to the 
whole atomic cell. This axis is identified with the spatial 
average of the magnetization density m n (f) in each cell 
n. This defines the local spin frame of reference. In this 
way we neglect the variation of the spin quantization axis 
within the cell during self consistency, avoiding the time- 
consuming numerical solution of the potential of coupled 
Schrodinger equations of the two spin channels. Within 
the local density approximation of density-functional the- 
ory, the exchange correlation potential has the same ref- 
erence frame as the local magnetization m„(r). Then for 
each atom we have a potential which is collinear in the 
local frame, and the solutions of the Schrodinger equa- 
tion, R l °l s (f;E) and H l °l s (r;E), depending on the spin 
index s of the local frame. 

The solution of the Schrodinger equation separately for 
each spin channel provides also the diagonal i-matrix of 
each atomic cell n in the local frame of reference: 



4T(£) 











(8) 



Then the i-matrix is rotated from the local to the global 
spin frame of reference using the spin rotation matrix 



tt\E) = U n t^(E)Ul 



U n being given by 



U r . 



cos(%0e-5' 
. sin(%-)e^ 



-sin(^)e 
cos(%0e^ 



(9) 



(10) 



The polar angles 8 n and <p n define the direction of the 
local magnetic moment with respect to the global spin 
frame of reference. Normally, 9 n and <j) n vary within 
the atomic cell, but in the approximation used here, av- 
erage angles are defined for each cell via an averaging 
of the magnetization density within the cell. Of course, 
when self-consistency is achieved, both the averaged and 
the point-by-point varying magnetization direction can 
be extracted from the output density matrix. Thus the 
assumption of a unique spin direction in each cell is only 
made for the spin-dependent potential. 

The i-matrix in the global spin frame of reference can 
be rewritten in the following way: 



J 



tf>\E) = U n 



-(tfr(25) + t^(E))l + -(t^(E) - t^(E))a z 



Ut 



(11) 



r 



with <t z is the z component of the Pauli matrices: 



1 

-1 



It is convenient to define the projection matrices <r ns 
for the local spin-up (f) and spin-down (J.) directions as: 



(12) 
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= 2^ 



(l±(T z )Ul = {cr ns y (+ for s 



Then t^ oh (E) is written as 

tt h {E)=t%{E)a 



=t, - for a =[) 
(13) 



(14) 
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In the collinear case the local and global frames are iden- 
tical and the projection operators reduce to: 



(collinear case). (15) 



At this stage, the difference between the ^-matrices 

At^ oh = £S lob _ t n is calculated in order to get all the 
ingredients to solve the Dyson equation for the structural 

o 

Green function (t„ has been defined in the global frame 
in Eq. ©). This is the analogue of Eq. (J2J) in matrix 
form in spin space: 

G stI (E) = G stT (E) + G stI (E)At sloh (E)G str (E). (16) 



Here, in analogy to eq. and eq. ©, G s t r (E) are matri- 
ces of size 2 x 2 in spin space, size (l max + 1) 2 x (l m ax + 1) 2 
in angular momentum space, and size N x N (with N the 
number of sites) in real space; all these indices are com- 
bined to form 2 x (Z max + l) 2 X iV-dimensional matrices. 
The i-matrix itself is diagonal in real space site indexes. 
The solution of Eq. (f ICift for the structural Green function 
requires matrix inversion, yielding G s t r (E) in the global 
frame: 



G stI (E) = G stI (E)(l - At^ ob (E) G stI (E))~\ (17) 



Equation Q can be now rewritten in the non-collinear 
case in order to obtain the Green function in the global 
frame. Using the matrices cr ns (Eq. |Q ) to project the 
local wavefunctions to the global frame, the Green func- 
tion is written as: 



G glob (i?„ + r, R n , + E) = -iV^ifefe; £)ff&.(f>; E)a ns 

Ls 

+ E Rtlsi^E^Gi^iE^s-R^A^^E). (18) 

LL'ss' 
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If needed, the Green function can be rotated to the local 
frame of any atom by the use of the transformation ma- 
trices U n (Eq. ^lOjl l. We point out that, even in the local 
frame of reference, the Green function is not in general 
diagonal in spin space. Finally we calculate the charge 
density and spin density from Equations (|3I4II . The spin 
dependent local density of states within the Wigner-Seitz 
cell WS of each site n is: 

n ns (E) = -- / lmG ss {R n + f, R n + r*; E)d 3 r (19) 
K Jws 

The spin density m = (m x ,m v ,m z ) (Eq.Q) is non- 
collinear. The new polar angles at each site n can then 
be obtained for each point by 

tan0„(f) = tan0„(r) = ^|t5 ( 2 0) 

m n (r) mn\ r ) 

or as an average over the local Wigner-Seitz cell 

fwo m « (r)df f WQ m% (r)df , x 

tan(9„ = J „ ws tan0 n = r WS Z,i ( 21 ) 



/ws m n(r)dr' 



Iws m n(r)df' 



III. APPLICATIONS 

As an application of our method, we study the mag- 
netic state of 3d-atom clusters in and on the Ni(001) 
surface. In a first step, we study the adatom proper- 
ties, which are already known from previous work. In a 
second step, we perform calculations for 3d dimers and 
trimers and use the understanding gained from the single 
adatoms in order to explain the results. 

Our calculations henceforth are based on the Local 
Spin Density Approximation (LSDA) of density func- 
tional theory with the parametrization of Vosko et alM^. 
The full nonspherical potential was used, taking into ac- 
count the correct description of the Wigner-Seitz atomic 
cells^ 3 . Angular momenta up to Z max = 3 were included in 
the expansion of the Green functions and up to 2^ max = 6 
in the charge density expansion. Relativistic effects were 
described in the scalar relativistic approximation. 

First, the surface Green functions arc determined by 
the screened KKR method 2 ^ for the (001) surface of Ni 
which serves as the reference system. The equilibrium 
lattice parameter of Ni was used (6.46 a.u. « 3.42 A). 
To describe the impurities on the surface (later we re- 
fer to these as adatoms and to the impurities sitting in 
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the first surface layer as inatoms), we consider a clus- 
ter of perturbed potentials which includes the potentials 
of the impurities and the perturbed potentials of several 
neighboring shells, with typical size ranging from 19 per- 
turbed sites for the single impurity to 32 for the dimers 
and trimers; in all cases, at least the first neighboring 
sites of the impurity atoms were taken into account in 
the calculation to ensure the correct screening of the im- 
purity potential. Test calculations have shown that this 
is adequate for our work; this is a merit of the Green func- 
tion method, in which the correct boundary conditions of 
the host (in our case of the host surface) are included in 
the Green function via the Dyson equation. We consider 
the adatoms at the unrelaxed hollow position in the first 
vacuum layer, and the inatoms at the unrelaxed position 
in the first surface layer. 

The orientations assigned to the spin moments of the 
impurities are always relative to the orientation of the 
substrate moment, which we take as the global frame. 
This, in turn, depends on delicate physical quantities 
such as the magnetic anisotropy energy, which cannot be 
related to the local properties of the small clusters that 
we study. In the present approach such effects arising 
from spin-orbit interaction are not included. The direc- 
tion of the host moments must therefore be considered as 
an input parameter from experiments or from indepen- 
dent ab-initio calculations. 



A. 3d single adatoms and inatoms 

3d adatoms on Fe(OOl) and on Ni(OOl) have 
been already studied previously, using the KKR 
me |-j 10(: j25 i 26 4 27 J 28 j n |_ ne atomic sphere approximation. 

Here we repeat the calculations of 3d adatoms on Ni(OOl) 
using the full potential method (a detailed work on Fe 
and Co on Ni(OOl) is presented in a recent article^). We 
give a brief analysis of the results, which are basically un- 
changed, in order to use them as a step for understanding 
the behavior of dimers and trimers later on. 

A collinear calculation of the magnetic state of a single 
adatom on a ferromagnetic substrate can give in some 
cases two solutions: one with ferromagnetic coupling 
(FM) to the substrate and one with antiferromagnetic 
coupling (AF). One of these states will correspond to the 
real ground state, and the other to a local minimum; this 
is actually a local minimum with respect to collinear vari- 
ations of the magnetic moment, since the angle between 
the local moment and the substrate moment cannot be 
varied in a collinear calculation. From total energy cal- 
culations of the two states, the ground state can be then 
determined. In some cases, when the intra-atomic ex- 
change field is not strong (beginning or end of the 3d 
series), only one of the two minima exists. On the other 
hand, if non-collinear effects are included in the calcula- 
tions, one of the two minima usually becomes unstable 
against an angular rotation of the moment, i.e., it is then 
actually a saddle point. 
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FIG. 1: 3d adatoms and inatoms on Ni(OOl): (a) Energy dif- 
ference between the AF and FM coupling, the values related 
to adatom and inatoms are described by respectively full and 
empty black diamonds; (b) magnetic moments of the adatoms 
(black triangles) and inatoms (green circles) within the 2 pos- 
sible magnetic configurations FM and AF; (c) the variation 
of the magnetic moments of Ni first Nearest Neighbors of the 
adatoms. 



The full diamonds in the Fig. Ufa) show the energy 
difference between the AF and the FM solution for 3d 
adatoms on Ni(OOl). The first elements of the 3d se- 
ries (Sc, Ti, V, Cr) are AF coupled to the substrate 
whereas the coupling of Mn, Fe, Co and Ni is FM. Sc 
(AF), Ti (AF), and Ni (FM) are characterized by a sin- 
gle solution. Clearly, the AF-FM transition occurs when 
the adatom atomic number changes from Cr [Z = 24) 
to Mn (Z = 25). This transition can be interpreted as 
in the case of the interatomic interaction of magnetic 
dimers in terms of the energy gain due to the for- 
mation of hybrid states with the Ni substrate as the 3d 
virtual bound state comes lower in energy with increasing 
Z. An explanation (see FigEJ can be given in terms of 
the d-d hybridization between the adatom 3d states and 
the Ni substrate 3d states. Energy is gained when a half- 
occupied d virtual bound state (VBS) at Ep is broadened 
by hybridization with the Ni minority 3d states, which 
lie at Ef (the Ni majority d states are fully occupied 
and positioned below Ef). For the early 3d adatoms 
(Fig(2Ja), it is the majority d VBS which is at Ep, thus 
the majority-spin direction of the adatom is favourably 
aligned with the minority-spin direction of Ni, and an 
AF coupling arises. For the late 3d adatoms (Fig|2t), 
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on the contrary, the minority d VBS is at Ep, and this 
aligns with the Ni minority d states; then a FM coupling 
arises. For our purposes we keep in mind that, since Cr 
and Mn are in the intermediate region, i.e., near the AF- 
FM transition point, their magnetic coupling to the Ni 
substrate is weak; this has consequences to be seen in the 
behavior of dimers, trimers, etc., in the next subsections. 

The magnetic moments of the adatoms and Ni first 
neighbors in the surface layer are shown in Fig. ^b) 
and (c). Evidently the moment of the Ni first neigh- 
bors is strongly affected by the adatoms. Especially in 
the AF state for Mn, Fe, and Co adatoms, the Ni mo- 
ment is strongly reduced and the FM configuration is 
stable. As regards the adatom moments, due to its half 
filled d band the Mn adatom carries the highest mag- 
netic moment (4.09 (xb) followed by Cr (3.48 /is) and Fe 
(3.24/xb). 

To understand the effect of coordination and stronger 
hybridization on the magnetic behavior of the adatoms, 
we take the case of impurities sitting in the first surface 
layer (inatoms). We carried out the calculations for V, 
Cr, Mn, Fe and Co impurities. The corresponding spin 
moments are shown in Fig^D (green circles), and the 
FM-AF energy differences are shown in Fig. (open 
diamonds and dashed line). 

Compared to the adatom case, the spin moments are 
reduced, especially for V and Cr. This effect is expected 
due to the increase of the coordination number from 4 
to 8 and the subsequent stronger hybridization of the 3d 
levels with the host wavefunctions. Moreover, the energy 
difference AE between the AF and FM solutions is af- 
fected. The trend can be understood as follows. In the 
case of Cr, the reduction of the local magnetic moment M 
is accompanied by a reduction of the exchange splitting 
AEx as AEx ~ I'M, where / w leV is the intra-atomic 
exchange integral. This means that, for the inatom, the 
occupied 3d states are closer to Ep than for the adatom. 
In turn, this intensifies the hybridization of these states 
with the Ni 3d states (which are close to Ep). At the 
same time, also the higher coordination number inten- 
sifies the d-d hybridization. The hybridization-induced 
level shift in the AF configuration increases, and the en- 
ergy of the AF state is thus lowered. The same mecha- 
nism is responsible for the weakening of the FM coupling 
of Mn inatom compared to the adatom. Similarly, the 
stronger hybridization of the Co-inatom d-states stabi- 
lizes even more its FM configuration due to the energy 
gain from the broadening of the d virtual bound state. 



B. Adatom and inatom dimers 

Having established the single adatom behavior, we 
turn to adatom dimers. We considered three geometries 
of increasing distance: dimers as first, second, and fourth 
neighbors. We will discuss the magnetic interaction be- 
tween the dimer atoms and the resulting magnetic order, 
first looking only at collinear states and then allowing 



for non-collinear order. We will see how, in certain cases, 
the collinear state reduces the symmetry, while the non- 
collinear state restores the full symmetry of the system. 
Non-collinear order is finally established for certain first- 
neighbor dimers. 

Fig.|3|represents schematically the different considered 
geometrical configurations of impurity dimers residing 
on the surface. We have investigated the dimer-l-type 
of geometry (the adatoms are first neighboring atoms), 
dimer-2-type (the adatoms are second neighbors) and 
dimer-3-type (the adatoms are fourth neighbors). This 
allows us to monitor the strength of the magnetic cou- 
pling as a function of the distance. Three collinear mag- 
netic configurations were treated: (i) antifcrromagnetic 
coupling within the dimer leading to a ferrimagnetic solu- 
tion (Ferri), (ii) ferromagnetic coupling within the dimer 
with both atoms ferromagnetically coupled to the sub- 
strate (FM), or (iii) ferromagnetic coupling within the 
dimer with both atoms antiferromagnetically coupled to 
the substrate (AF). 

Our calculations include V, Cr, Mn and Fe dimers. 
We found that all V and Fe dimer types behave like 
the adatoms: in all geometries, both V atoms are AF 
and both Fe atoms are FM. On the other hand, Cr 
and Mn dimers show magnetic frustration. As shown 
in Fig. both the Cr-dimer-1 and Mn-dimer-1 show 
(in a collinear calculation) a Ferri ground state (see Ta- 
ble HJ. With increasing distance between the adatoms, a 
transition occurs to the single adatom magnetic behavior 
which is AF for Cr-dimers and FM for Mn-dimers. It is 
clear that, in the dimer- 1 case, there is a competition of 
exchange interactions. 

When we allow for a rotation of the magnetic moments, 
non-collinear solutions are obtained for the Cr- and Mn- 
dimer-1 systems. On the other hand the magnetic cou- 
pling of the V- and Fe-dimer-1 remains collinear. Let us 
start with Cr-dimer-1: Fig.^Ja) represents the collinear 
magnetic ground state. As one expects from the adatom 
picture, both adatoms forming the dimer tend to cou- 
ple AF to the substrate but due to their half filled d 
band they also tend to couple AF to each other. This 
can be understood in terms of the Alexander-Anderson 
modeli22*2i To give a short explanation (see Fig. [20, both 
Cr and Mn have their majority-spin VBS occupied, be- 
low Ep, and the minority-spin VBS unoccupied, above 
Ep. An antiparallel configuration between the moments 
in a Mn or Cr dimer lowers the energy, because the occu- 
pied d VBS of each atom hybridizes with the unoccupied 
d VBS of the other atom and is shifted to lower eigenval- 
ues. Contrary to this, a parallel coupling does not lower 
the energy, since there is no level shifting, but only level 
broadening of the majority d-VBS. Since these are fully 
occupied, the broadening brings no energy gain. 

Thus there is a competition between the interatomic 
coupling within the dimer, which drives it to a Ferri state, 
and the exchange interaction with the substrate, which 
drives the moments of both atoms in the same direc- 
tion: AF for Cr and FM for Mn. As discussed in the 
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V 
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DOS 
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J broadening 



b. Fe-Ni: FM 
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shift A 



hybridi-iy 
zation /V 



shift A 



FIG. 2: Alexander-Anderson model for neighboring magnetic atoms: (a) Early 3d transition elements in interaction with Ni 
surface atoms; (b) Late 3d transition elements in interaction with Ni surface atoms; (c) Cr or Mn dimer. 





Cr 


Mn 




Dimer 1 Dimer 2 Dimer 3 


Dimer 1 Dimer 2 Dimer 3 


Efm — -E'p e rri(eV) 
E\F - -BFcrri (eV) 


0.451 0.130 0.120 
0.433 -0.093 -0.112 


0.065 -0.242 -0.239 
0.496 0.187 0.233 



TABLE I: Energy differences between the Ferri solution and the FM (AF) configuration for the three types of dimers investigated. 



previous subsection, the magnetic exchange interaction 
(MEI) to the substrate is relatively weak for Cr and Mn. 
Thus, the intra-dimer MEI is stronger than the MEI 
with the substrate, and in the collinear approximation 
the ground state is found ferrimagnctic (Ferri). Remov- 
ing the collinear constraint, a compromise can be found 
such that both atoms are AF coupled to each other and 
at the same time (for Cr) slightly AF coupled to the 
substrate. This is shown in Fig^Jb): the Cr adatom 
moments are aligned antiparallcl to each other and ba- 
sically perpendicular to the substrate moments. How- 
ever, the weak AF interaction with the substrate causes 
a slight tilting towards the substrate, leading to an angle 
of 94.2° instead of 90°. We also observe a very small tilt- 
ing (« 0.3°) of the magnetic moments of the four outer 
Ni atoms neighboring the Cr dimer (the two inner Ni 
atoms do not tilt for symmetry reasons). 

Despite the above considerations, the collinear Ferri 
state (Fig^a)) is also a self-consistent solution of the 
Kohn-Sham equations, even if the collinear constraint 
is removed. Total energy calculations are needed in or- 
der to determine if the non-collinear state is the true 
ground state, or if it represents a local minimum of energy 
with the collinear result representing the true ground 
state. After performing such calculations we find that 
the ground state is collinear with an energy difference 



of Ai?Ncoi-Forri = 39.84 meV (increasing the angular- 
momentum cutoff to i max = 4 brought no significant 
change to this result). 

The case is different for Mn dimers. Fig. [5] shows 
the collinear and the non-collinear solutions. The 
dimer atoms couple strongly antiferromagnetically to 
each other but the single Mn adatoms tend to cou- 
ple (weakly) ferromagnetically to the substrate. Both 
adatom moments, while aligned AF with respect to each 
other, are tilted in the direction of the substrate mag- 
netization, as opposed to the Cr-dimer. With a ro- 
tation angle of ss 72.6°, the tilting from the 90° con- 
figuration is rather large. Also the Ni moments are 
tilted by 7.4°. The main difference with the case of Cr- 
dimer-1 is that for Mn-dimer-1 the non-collinear solu- 
tion is the ground state (total energy calculations yield 
Ai^Ncoi-FciTi = —13.45 meV). The spin moments of the 
V, Cr, Mn, and Fe dimers are given in Table UTTl 

In both cases (Cr and Mn dimers) the frustrated 
collinear solution is asymmetric, while the non-collincar 
ground state restores the twofold symmetry of the sys- 
tem. The differences in energy between the Ferri and the 
non-collinear solutions are small and can be altered ei- 
ther by using a different type of exchange and correlation 
functional such as GGA or LSDA+C7, or after relaxing 
the atoms. We note, however, that in a test calculation 



DimeM Dinner 2 




Mn dimers Fe dimers 

onNi(001) onNi(001) 

FIG. 3: Different geometrical configurations considered for 
dimers at the surface of Ni(OOl). Dimer-l-type corresponds to 
the case where the atoms are first neighboring atoms, dimer- 
2-type where the atoms are 2'NN and finally dimer-3-type to 
4'NN. The collinear magnetic ground state are also shown for 
V, Cr, Mn and Fe dimers. 




FIG. 4: Most stable configurations of Cr dimer-l-type ob- 
tained with (a) the collinear KKR method and (b) the non- 
collinear KKR method. The rotation angle with respect to 
the z axis is equal to 94.2°. The collinear state is the ground 
state, with the non-collinear state being a local minimum (see 
text). 



we found the Cr single-adatom relaxation to be small 
(3.23 % inward with respect to the interlayer distance), 
and thus we believe that the relaxation cannot affect the 
exchange interaction considerably. 

As a cross-check, it is interesting to compare these non- 
collinear ab-initio results to model calculations based on 
the Heisenberg model with the exchange parameters fit- 
ted to the total energy results. We assume a classical 
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FIG. 5: Most stable configurations of Mn dimer-l-type ob- 
tained with the collinear KKR method (a) and non-collinear 
KKR method (b). The rotation angle with respect to the z- 
axis is equal to 72.6°. The non-collinear state is the ground 
state. 



spin Hamiltonian of the form 

H = -\Y t Ji&e j . (22) 

Here, e is a unit vector defining the direction of the mag- 
netic moment and i and j indicate the dimer atoms and 
their first Ni neighbors. We can evaluate the interatomic 
exchange constants Jcr-Ni, ^Mn-Ni, Jkn-Mn and Jcr-Cr 
via a fit to the total energy obtained from collinear LSDA 
calculations of the FM, AF, and Ferri configurations. 
Taking into account only first-neighbor interactions and 
neglecting the rotation of Ni moments, we rewrite the 
Hamiltonian for the dimer in terms of the tilting angles 
6*i and 9-2 of the two Cr (or Mn) atoms (the azimuthal 
angles <fr do not enter the expression because of symmetry 
reasons): 

H = — Jcr-Cr COs($i — #2) — 4 J(>- Ni(cOS 01+COS #2)+C0nst. 

(23) 

We note the two extreme cases arising from this Heisen- 
berg Hamiltonian: (i) |Jc r -Ni| ^> |Jcr-Cr| leads to 
the stabilization of the collinear FM or AF configura- 
tion (adatom-like behavior) and (ii) | Jo- Nil | Jcr-Cr| 
leads to antiferromagnetic coupling within the dimer if 
Jcr-Cr < 0. Within the Heisenberg model the Ferri so- 
lution and the non-collinear solution with 8 = 90° have 
the same energy. 

Table [H] summarizes the estimated exchange param- 
eters. Two effects are striking: (i) The strong an- 
tiferromagnetic Cr-Cr and Mn-Mn interaction for the 
dimer-1 (nearest neighbors), being an order of magni- 
tude larger than the exchange interactions with the sub- 
strate and being responsible for the stabilization of the 
non-collinear state structures shown in Figs.0]and[5] (ii) 
The very weak Cr-Cr and Mn-Mn interactions in the 
dimer-2 and -3 configurations. Whereas for the nearest- 
neighbors configuration (dimer-1) the direct overlap of 
the rf-wavefunctions of the Cr and Mn atoms leads to 
the strong coupling, this overlap is missing for larger dis- 
tances and the interaction can only proceed through the 
substrate. However, this interaction is weak, in fact con- 
siderably smaller than the interaction of both adatoms 
with the four neighboring Ni atoms of the substrate. 
Therefore these dimers are effectively decoupled, and be- 
have like the isolated adatoms, being antiferromagneti- 
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(a) 


00 


Jij (meV) 


Dimer 1 


Dimer 1 


Dimer 2 Dimer 3 


JCr-Ni 


-1.3 


-11.6 


-13.9 


-14.5 


Jcr-o 


-189.1 


-221.3 


-9.2 


-2.0 


^Mn-Ni 


13.0 


27.0 


26.8 


29.5 


^Mil-Mil 


-138.2 


-140.2 


13.7 


1.5 



TABLE II: Values of magnetic exchange parameters Jij for 
Cr and Mn dimers on Ni(001), fitted from collinear first- 
principles total energy calculations (b) and obtained by the 
Lichtenstein formulae- (a) (Jcr-Ni and Juin-Ni are averaged 
over the different Ni first neighbours of the dimer atoms). 
Positive Jij correspond to ferromagnetic interactions, nega- 
tive Jij to antiferromagnetic ones. 



ent kind of coupling that each inatom has with the the 
substrate (AF or FM). One notices also that the mag- 
netic moments in the ground state decrease compared 
to the values obtained for the single inatoms and sin- 
gle adatoms. When the rotation of the moments is al- 
lowed, Cr-dimer can be stabilized at an angle of 107° 
(instead of 94.2° found for the adatom-dimer case), and 
Mn-dimer at an angle of 80.9° (instead of 72.6°). Thus 
the non-collinear solutions obtained for inatom-dimers 
are rather similar to what was obtained for adatom- 
dimers. Energetically, however, both the Cr and the Mn 
inatom dimers show a lower total energy in the collinear 
Ferri state (for Cr, AI?Ncol-Ferri = 24.11 meV; for Mn, 
A-ENcol-Ferri = 22.5 meV). 



cally coupled to the substrate in the case of Cr and fer- 
romagnetically for Mn. The exchange constants Jy fit- 
ted to total energy results can be compared to the ones 
obtained by using the Lichtenstein formula^ (starting 
from the Ferri ground state). This rests on the force 
theorem, and yields the exchange constants relevant to 
an infinitesimal rotation of the moments. The results of 
the two methods agree best for the Mn-Mn interaction, 
and reasonably well for the Cr-Cr interaction, but not for 
Mn-Ni and Cr-Ni. 

With the parameters from Table [H] one can also re- 
calculate the non-collinear structure of the ground state. 
The agreement with the ab-initio results is quite reason- 
able. For the Cr dimer, one finds a slightly smaller tilting, 
i.e. 96° instead of 94.2°, while for the Mn dimer the angle 
is 67.3° instead of 70.6°. 

The differences in energy calculated within this sim- 
ple model, show that the Cr-dimer-1 has a non-collinear 
ground state (Ai/Ncoi-Ferri = —9.7 meV) as well as the 
Mn-dimer- 1 (Ai?Ncoi- Ferri = —41.6 meV). The discrep- 
ancy obtained for the case of Cr-dimer-1 (the LSD A cal- 
culation gives the collinear Ferri ground state) can be at- 
tributed to the restrictions of the Heisenberg model. For 
instance, for the Ferri and non-collinear configurations, 
the Cr moments are slightly different, and also the re- 
duction of the Ni moments as a function of the rotation 
angle {e.g. for the single adatom) cannot be described 
by the Heisenberg model, where the absolute values of 
the moments are assumed to be constant. Within the 
Heisenberg model, the Ferri solution (with Q\ = 0° and 
6-2 = 180°) is degenerate with the non-collinear solution 
(#i,2 = 90° with AF coupling within the dimer). 

To evaluate the effect of change in coordination and 
hybridization, we have undertaken a study of inatom 
first-neighbor dimers for V, Cr, Mn and Fe. The V and 
Fe inatom-dimers were found to behave like the adatom 
dimers. The V dimer is in an AF state, the Fe dimer 
in a FM state, while the Cr and Mn dimers arc in a 
Ferri state (in case of collinear constraint). The spin mo- 
ments in the collinear and non-collinear states are given 
in Table UTTl Within the Ferri-dimers, the difference be- 
tween the moments of the two atoms is due to the differ- 



C. Trimers 




FIG. 6: Non-collinear state of the Mn trimer on Ni(001) sur- 
face. Side view (a) and front view (b) are shown. This rep- 
resents a local minimum in energy, with the collinear state 
being the ground state (see text). 



Following the same procedure as for dimers we first 
investigated several collinear magnetic configurations for 
the most compact trimer on the Ni(001) surface, which 
has the shape of an isosceles rectangular triangle (see 
Fig. 0) of side y/2a/2 and hypotenuse a (a is the Ni fee 
lattice constant). It is expected, and verified by total- 
energy calculations, to find the jfj, configuration as the 
collinear magnetic ground state for Cr and the | J,f for the 
Mn trimer (f means an atomic moment parallel to the 
substrate, J, an antiparallel one; the middle arrow repre- 
sents the direction of the atomic moment at the right- 
angle corner of the triangle). In Table ITVl the energy dif- 
ferences among the possible collinear configurations are 
given; for the Iff and J, J, J, Cr trimers, no self-consistent 
solution could be found. 

Allowing free rotation of the magnetic moments leads 
to no change for the Cr trimer ||| — the state remains 
collinear (within numerical accuracy). On the other 
hand, for the Mn trimer a non-collinear solution is found 
(Fig.|SJ) with the nearest neighbours almost antiferromag- 
netic to each other, but with a collective tilting angle 
with respect to the substrate. This tilting angle is in- 
duced by the ferromagnetic MEI between the central Mn 
atom with the substrate, competing with the antiferro- 
magnetic MEI with its two companions. The top view 
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Dimertype V(AF) Cr(Fcrri) Cr(Ncol) Mn(Fcrri) Mn(Ncol) Fc(FM) 

OnNi(OOl) (-1.28,-1.28) (-3.04,3.05) (3.03,3.03) (-3.84,3.69) (3.75,3.75) (3.10,3.10) 
InNi(OOl) (-0.32,-0.32) (-2.00,1.96) (1.97,1.97) (-3.32,3.20) (3.26,3.26) (2.88,2.88) 

TABLE III: Atomic spin moments (in /is) of the adatom and inatom dimers (of type 1, i.e., nearest-neighbors) in the collinear 
and non-collinear configurations. A minus sign of the collinear moments indicates an antiparallel orientation with respect to 
the substrate magnetization. Embedding the dimer into the surface causes, as expected, a decrease of the spin moments due 
to stronger hybridization of the d wavefunctions. 



Magn. Config. m ITI UT TIT TTT III 
Cr: E-E ln {eV) 0.420 0.390 0.193 
Mn: E-E ln {eV) 0.116 0.318 -0.184 0.239 0.817 

TABLE IV: Energy differences between the different calcu- 
lated collinear magnetic configurations with the ||| configu- 
ration. The direction of the arrows represents the direction of 
the atomic moments relative to the substrate magnetization 
(t parallel, [ antiparallel). The middle arrow represents the 
atom at the right-angle corner of the trimer. 

of the surface shows that the in-plane components of the 
magnetic moments are collinear. 

The tilting is somewhat smaller (21.7°) for the two 
Mn atoms with moments up than for the Mn atom 
with moment down (28.5°). Also the neighboring Ni- 
surface atoms experience small tilting, with varying an- 
gles around 4° — 10°. From the energy point of view, the 
ground state is the collinear one, |||, with an energy dif- 
ference of Ai?Ncoi-TiT = 22.92 mcV with respect to the 
non-collinear solution. 

We have also investigated the cases where the trimers 
are sitting in the surface layer. No non-collinear solu- 
tion was found, while there is no change in the collinear 
ground state which is ||| for the Cr trimer and TIT for 
Mn trimer. 

One should note that the moments of the two first 
neighboring impurities are almost compensated in the 
Fcrri solution. The third moment determines the total 
interaction between the substrate and the trimer which 
has then a net moment coming mainly from the addi- 
tional impurity. This interaction is identical to the single 
adatom (or inatom) type of coupling. 

IV. SUMMARY 

We have presented a formalism for the treatment of 
non-collinear magnetic clusters on surfaces and in bulk, 
based on the Green function technique of Korringa, Kohn 
and Rostoker, and on spin density functional theory. We 
have applied the formalism on the study of small tran- 
sition metal clusters (dimers and trimers) on and in the 
Ni(001) surface. 

Emphasis was placed on Cr and Mn clusters, for which 
we found that magnetic frustration can lead to non- 



collinear magnetic order. The origin of the frustration is 
the competition of the antiferromagnetic exchange cou- 
pling among the Cr or Mn atoms with the antiferromag- 
netic (for Cr) or ferromagnetic (for Mn) exchange cou- 
pling between the adatoms and the substrate. In this 
respect, the result is different than the prototype non- 
collinear configurations arising from antiferromagnetic 
interactions among atoms in triangular geometry. 

We found that Cr and Mn first-neighbouring adatom 
dimers can show non-collinear behavior, while increasing 
the distance between the adatoms of the dimer leads to 
the same state as for single adatoms. The energy dif- 
ferences between the collinear ferrimagnetic state and 
the non-collincar one are A£^ ol _ Forri = 39.84 meV (the 

ground state is collinear), Ai^coi-Fcrri = ^13-45 meV 
(the ground state is non-collinear). Embedding the 
dimers in the first surface layer restores the Ferri collinear 
solution as a ground state also for Mn adatom dimers ( 
A^oi-Ferri = 24.11 meV, A<« ol _ Ferri = 22.5 meV). 

Our ab-initio results for dimers are compared to the 
solution of a classical Heisenberg model with exchange 
parameters fitted to total energy calculations. The agree- 
ment for the tilting angles in the non-collincar state is 
good, but the Heisenberg model does not capture the 
collinear ground state for the Cr dimer. This discrep- 
ancy occurs because the Heisenberg model is restricted 
to constant absolute values of localized spins. 

The trimers studied so far are characterized by a 
collinear ground state: ||| for the Cr trimer and TIT 
for the Mn trimer. The Mn trimer has also a non- 
collincar metastable solution with an energy difference 
Ai?Ncoi-TlT = 22.92 meV. 

We believe that the energetic proximity of the collinear 
to the non-collinear states is directly related to the weak- 
ness of the exchange interaction with the Ni substrate. 
Replacing it by an fee Fe substrate will possibly change 
the ground state drastically. Work in this direction 
progress and will be reported elsewhere. 
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